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Abstract 



We consider the problem of optimal investment and consumption in a class of multidi- 
mensional jump-diffusion models in which asset prices are subject to mutually exciting 
jump processes. This captures a type of contagion where each downward jump in an 
asset's price results in increased likelihood of further jumps, both in that asset and in the 
other assets. We solve in closed-form the dynamic consumption-investment problem of 
a log-utility investor in such a contagion model, prove a theorem verifying its optimal- 
ity and discuss features of the solution, including flight-to-quality. The exponential and 
power utility investors are also considered: in these cases, the optimal strategy can be 
characterized as a distortion of the strategy of a corresponding non-contagion investor. 
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1 Introduction 



The recent financial crisis has emphasized the relevance of jumps for understanding the various 
forms of risk inherent in asset returns and their implications for asset allocation and diver- 
sification. Most portfolios, from those of individual investors to those of more sophisticated 
institutional investors including University endowments, suffered badly during the latest crisis 
episode, with many commonly employed asset allocation strategies resulting in large losses in 
2008-09. Reasonably diversified portfolios can survive a single isolated negative jump in asset 
returns. However, jumps that tend to affect most or all asset classes together are difficult to 
hedge by diversification alone. Moreover, additional jumps of this nature seemed to happen in 
close succession, as if the very occurrence of a jump substantially increased the likelihood of 
future jumps. 

Motivated by these events, we consider in this paper the issue of optimal portfolio con- 
struction when assets are subject to jumps that share the qualitative features experienced most 
vividly during the recent financial crisis. These salient features include the fact that multiple 
jumps were observed, at a rate that was markedly higher than the long term unconditional 
arrival rate; these jumps affected multiple asset classes and markets; and they affected them 
not necessarily at the same time, but typically in close succession over days or weeks. 

To capture these key elements, we consider a model for asset returns where a jump in 
one asset class or region raises the probability of future jumps in both the same asset class 
or region, and th e other classes or regions. Jump processe s of this type were first intro- 



duced by_ 



Hawked (Il97ll ) with further developments due to iHawkes and Oaked (Il974[ ) and 



OakesI fll975h . Mode 



seismology (see, e.g., 



s of this t ype 



lave been used in epidemiolog y, neurophysiology and 



BriUingeil (119881) and lOeata and Akaikel (11982 )' 



Reynaud-Bouret and SchbathI (120 lOl )). credit deri vatives 



action times and price change s at high frequency (IBowshei 



jrais et al 



, gen ome analysis (see 



(120 lOl )). to model trans- 



(l2007l )). trading activity at different 



maturities of the yield curve (ISalnaon and ThamI ( 120071 )) and propagation phenomena in social 
interactions (ICrane and Sornettd (120081 ).) 

We extend the pure jump Hawkes model employed in the above applications, in order to 
better represent financial asset returns. We add a drift to capture the assets' expected returns 
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and a standard, Brownian-driven, volatility component to capture their day-to-day normal 
variations. We call this m odel a "Hawkes jump- diffusion" by analogy with the Poisson jump- 



diffusion of iMertonI (Il976[ ). Unlike models typically employed in finance, the jump part of this 
model is no longer a Levy process since excitation introduces a departure from independence 
of the increments of the jumps. 

In the model, jump intensities are stochastic and react to recent jumps: a jump increases 
the rate of incidence (or intensity) of future jumps; running counter to this is mean reversion, 
which pulls the jump intensities back down in the absence of further excitation. In the univariate 
case, only "self excitation" can take place, whereas in the multivariate case "mutual excitation" 
consisting of both self- and cross-excitation (from one asset to another) can take place. We 
now illustrate the presence of the mutual excitation phenomenon by filtering jump intensities 
for jumps in the US financial sector stock index during the recent crisis. Figure [T] plots the 
estimated intensity of US jumps over time, filtered from the observed returns on indices of 
financial stocks in the US, UK, Eurozone and Asia. The excitation mechanism is apparent in 
the Fall of 2008, when jump intensities increase rapidly in response to each jump, most of them 
originating in the US, and to a lesser extent in the Winter of 2009, which looks more like a slow 
train wreck. 

The bottom panel of the plot shows that the filtered intensities contain information that 
is different from other measures of market stress, VIX and the CDS rate on financial stocks. 
In particular, VIX is a measure of total quadratic variation and as a result captures the total 
risk of the assets instead of just their jump risk. So the same jumps which cause the jump 
intensity to increase in the middle panel also cause VIX to increase in the bottom panel, but 
VIX includes Brownian volatility, making it a much noisier measure of jump risk. 

The purpose of the paper is to solve for the optimal portfolio of an investor who faces this 
type of risk in his/her investment opportunity set. By considering a more realistic model for 
jumps, incorporating mutual excitation, we are able to study the optimal portfolio of an in- 
vestor in a realistic setting where a jump that occurs somewhere will increase the probability 
of further jumps in other asset classes or markets. The model generates jumps that will tend 
to be clustered (as a result of the time series self-excitation), systematic (as a result of the 
cross-sectional excitation), but neither exactly simultaneous nor certain, since the excitation 
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Figure 1: Time- varying Jump Intensities During the Financial Crisis. The top panel shows 
time series of stock indices for the financial sector in the four regions, 2007-2009. The middle 
panel shows the intensity of US jumps filtered from the model, based on each 3a and above 
event identified as a jump. Each jump leads to an increase in the jump intensity, followed by 
mean reversion until the next jump. The bottom panel shows the time series of two alternative 
measures of financial distress, the VIX index and Markit's CDS index for the banking sector in 
the US. 
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phenomenon merely raises the probabihty of future jump occurrence. By analogy with epi- 
demics, the probability of getting infected increases in a pandemy but does not typically reach 
one, and there is an incubation period which can range in the case of financial markets from 
hours to days, depending upon subsequent news arrival, and once established the pandemy 
does not go away immediately. Furthermore, the model is multivariate and the contagion can 
be asymmetric, with jumps occurring in one asset class or market having a greater excitation 
potential for the other sectors or regions than jumps that originate elsewhere: for instance, 
most financial crises that originate in or transit through the US tend to have greater ramifica- 
tions in the rest of the world than crises that originate outside the US. Poisson jumps, whose 
intensities are constant, are not able to reproduce these empirical features, and this motivates 
our inclusion of the more general class of Hawkes jumps in the model. 



lios when asset returns can iump (see, e.g., 
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1984 
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1990), 
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kawa (1990 
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Han and Rachev (200i 
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OrtobeUi et al. (' 


2003). 
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(2001 
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2003) 
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novel aspect in the present paper is the inclusion of Hawkes jumps in asset returns: such 
jumps share the dual characteristics of being systematic, meaning that they affect multiple 
assets or asset classes at the same time, and mutually exciting, meaning that they affect the 
rate at which future iu mps occur in each asset class. By contrast, in the earlier model of 
Ait-Sahalia et al.l ( l2009l ). assets were subject to random jumps which could affect one or more 
asset or asset classes, but when they occurred, they were simultaneous and every asset in that 
sector or region would jump. Such jumps were also "Poissonian" , in the sense that the arrival 
of jumps today did not influence the future arrival of jumps. 

We show that the optimal portfolio solution in the model can be obtained in full closed-form 
in the log utility case, and in quasi- closed-form in other cases. Importantly, we show that the 
optimal solution becomes time-varying, with the investor reacting to changes in the i ntensity 
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of the jumps. For a log-utility investor, the solution remains myopic, as in the classical 
( I1971I ) problem with log-utility, in the sense that the investor does not need to take into account 
the full dynamics of the state variables. The log investor in our model holds at each point in 



time the same portfolio as a log investor who believes that jump intensities are constant, but 
his/her optimal portfolio weight is now constantly changing to reflect the time- variation in 
jump intensities. One consequence of this result is that each time a market shock occurs, the 
investor perceives an increase in jump intensities, and sells some amount of each risky asset, 
and invests the proceeds in the riskless asset, a behavior we interpret as a "flight to quality." 

Formal computations also work for both power and exponential utility investors, although 
we have not proved the appropriate verification theorem for these utilities. Nevertheless, the 
resultant investment strategy under contagion can be interpreted in terms of the equivalent 
strategy under the "non- contagion" assumption that jump intensities are constant. We find 
that under the contagion conditions of the model, the investor will choose a portfolio that is 
optimal for a non-contagion investor who has a specific distorted value of the intensities, which 
we characterize. This distortion of intensities has the effect of magnifying the investment in the 
risky assets: the contagion investor will go "longer" when the non-contagion investor is long 
in the risky asset, and will go "shorter" when the non-contagion investor is short in the risky 
asset. 

The paper is organized as follows. Section |2] presents the model for asset returns. Section 
|3] introduces the optimal portfolio problem when jumps are mutually exciting in the general 
case. Section H] specializes the solution to the case of an investor with log-utility and derives 
the optimal portfolio and consumption policy in closed-form, including a complete verification 
theorem that supports this policy. Section [5] develops some interesting market specifications 
which exhibit an explicit closed-form log-optimal policy. Section |6] explores some of the prop- 
erties of such explicit asset allocation policies. The exponential and power utility investment 
problems are outlined in Section [71 In these problems, the solutions can be characterized as 
distorted versions of the non-contagion solutions. Section [8] concludes. 

2 Mutually Exciting Jumps 

In this paper, jumps in one asset class not only increase the probability of future jumps in 
that asset class (self excitation) but also in other asset classes (cross excitation). In a mutually 
exciting model, the intensity of a jump counting process ramps up in response to past jumps. 
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A mutually exciting process is a special case of path-dependent point process, whose intensity 
depends on the path of the underlying process. Mutually exciting counting processes, Ni^i 
{I = 1, m), form an m-vector Nt = [Ni,t, Nm,t]' such thaj] 

P [Ni,t+At - Ni,t = MJ't] = XiAi + o (At) , (2.1) 
F[Nu+At-Ni^t>0\J't]=o{At) 

independently for each /. In the standard model specification we adopt in this paper, the 
Hawkes intensity processes = [Ai^i, A^,*]' have the integrated form 

Az,t = e-"'*Az,o + (1 - e-"'*)Az,oo + V" T V"'^*"'^^^.',.. / = l,...,m (2.2) 

^-^ Jo 

For all l,j = l,...,m the parameters Xi^oo,dij > and Xifl,ai > are constants. These 
parameter restrictions ensure the positivity of the intensity processes with probability one. 
Differentiation of equation (12. 2p shows that the intensity in asset class / has dynamics given 

by 

Em 
dijdNj,t. (2.3) 
i=i 

In other words, a jump dNj^g, occurring at time s G [0, t) in asset class j = 1, . . . , m, raises each 
of the jump intensities A/^f , / = 1, . . . , m, by a constant amount dij. The Zth jump intensity then 
mean-reverts to level A;^oo at speed ai until the next jump occurs. Equation 12.31 also reveals 
that [N, A) is a 2m-dimensional Markov process, while (A) alone is an m-dimensional Markov 
process. 

As a result, our model generates clusters of jumps over time, and jumps can propagate at 
different speed and with different intensities in the different asset classes depending on where 
they originate and which path they take to reach a given asset class or market. Free parameters 
in the model control the extent to which the two forms of excitation take place, the relative 
strength of the contagion phenomenon in different directions and the speed with which the 
excitation takes place and then relaxes. 

The model produces both cross-asset class and time series excitation. In the univariate 
self-exciting case, a typical sample path of one component of A is illustrated in Figure [2l Each 
jump increases the jump intensity, followed by mean-reversion. 

^In this paper, we work in a filtered probability space (ft, T , {J-'t)t>o, P) that satisfies "the usual conditions." 

6 




Jump Times 

Figure 2: Sample path of a Hawkes intensity, Xi^f 

This model also produces cross- asset, or mutual, excitation. Jumps in asset class / that 
occurred u units of time into the past raise the intensity of jumps in asset class j by dj^e""^", 
while conversely jumps in asset class j raise the intensity of jumps in asset class / by dije~'^'^. 
Jumps in a given asset class i also raise the intensity of future jumps in the same asset class. 
Figure |3] illustrates this with two assets. At time Ti there is a jump in the first asset class value, 
5*1. This jump self-excites the jump intensity Ai. This increase in Ai raises the probability of 
observing another jump in Si at the future time T2. These jumps have a contagious effect on 5*2 
since a jump in Si cross-excites the jump intensity of 5*2. This, in turn, raises the probability 
of seeing a jump in 5*2 at time T3. Latter on, at time T4, the jump in S2 raises the probability 
of seeing a jump in Si at some future time T5. The degree to which self- and cross-excitation 
matter in the model, and their relative strengths, is controlled by the parameters in d and ex.. 

We note that each compensated process Ni^t — Jq M,sds is a local martingale. A number of 
important additional properties hold for this model: 



1. Markov Generator: The Markov generator of this process acting on differentiable 
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T1T2T3 T4 T5 



Figure 3: Cross- and self-excitation in a two asset-class world: Sample paths of the jumps, asset 
prices and jump intensities. 
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functions g : (Z_|_)™ x 



— i- M is given by 



[Ag]{n,\) = Y, 



1=1 



dg 

ai (A/,oo - -^/) ^ + -^i ioin + ei,\ + di) - g{n, A)) 

OAi 



where ei = [6ii, Smi]' and di = [du, dmi]'- If 



E 



/ [Ag] {N„ A,) - [ai {Xe,oo - M 
Jo 



dg{N,,Xs 
dXi 



ds 



< oo 



for all t and g is different iable in A, the Dynkin formula says that for each t < T : 



E[g{NT,XT)\J't]=g{Nt,Xt)+E 



(2.4) 



2. Stationarity Assumption: Let the process A^ satisfy (12. 3p with Aq € M^; ctj > 0; dij > 
and with F = {Tij)ij=i^,,,^rn, where Tij = ajSij — dij, a positive (hence invertible) matrix 
where 6ij is the Kronecker symbol. Then the intensities A are stationary processes with 
bounded moments. Using the Dynkin formula, one can show that under this assumption 
the first moments fi(t) = E [A/(t)] converge as t — t- oo to the non- negative values fi{oo) = 
'^^=i{Sij — aJ~^dij)\j^oo- Similar formulas can be derived for the large time limits of the 
higher m oment functions. U sing these bounds, the ergodic theorem for semimartingales 



[see e.g. 



Khasminskiil ( 1l960l )) then implies that for any measurable function K{X) that 
satisfies a bound |i^'(A)| < M(l + ||Af ) 



lim ^ [ K{Xt)dt = lim E[K{Xt)] = [ K{X)ij{dX) 

a.s. where fi{dX) is the invariant (infinite time) measure of A*. 
3. AfRne Structure: The joint characteristic function has the affine form 



(2.5) 



Eo,Ao[e 



iuNx+ivXi'] 



exp 



Em 



where the deterministic functions A, Bi satisfy the Riccati equations 
dA 



^^^aiXi^ooBf, A(0) = 0; 
^ = -«zi?.-^(e-+^^-^'^^'^^-l); B,iO) = v,; 1 = 1,.... 



m 



(2.6) 
(2.7) 
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3 Optimal Portfolio Selection When Jumps Are Mutu- 
ally Exciting 

We now solve Merton's problem: the investor maximizes the expected utility of consumption by 
investing in a set of n risky assets and a riskless asset over the infinite time horizon t e [0, oo). 
The innovation is that the risky assets are now subject to shocks generated by an m-dimensional 
Hawkes jump-diffusion process. 

3.1 Asset Return Dynamics 

The riskless asset with price So,t is assumed to earn a constant rate of interest r > 0. The n 
risky assets with prices St = [Si^t, • • • , Sn,t]' follow a semimartingale dynamics with asset shocks 
generated by an m-dimensional Hawkes process. Specifically, we assume 
dSo,t 



= rdt, (3.1) 

n m 

= (r + Ri) dt + J2 <^ijdWj,t + Ji,iZi,tdNi,t, ^ = 1, n (3.2) 
j=i 1=1 

Here iV^ = [Ni^t, .... Nm,t]' is an m— dimensional, m < n, vector of mutually exciting Hawkes 
processes with intensities = [Ai^^, . . . , Xjn,t]' that follow the Markovian dynamics 

Em 
. dijdNj^t, / = 1, . . . ,m, (3.3) 

with constant parameters ai > 0, A;,oo > 0, and dij > 0. Under the condition that F is a 
positive matrix, the A process is stationary. 

The vector Wt = [Wi,t, ■ ■ ■ , Wn,t\' is an n— dimensional standard Brownian motion. Ji^iZi^t 
is the response of asset i to the lih. shock where Zi^ti a scalar random variable with probabil- 
ity measure i'i{dz) on [0,1], is scaled on an asset-by- asset basis by the deterministic scaling 
factor Ji^i e [—1,0]. For clarity, we chose to include only negative asset jumps in the asset 
price dynamics, since those are the more relevant ones from both a portfolio risk management 
perspective and their contribution to mutual excitation. 

We assume that the individual Brownian motions, the Hawkes process and the random vari- 
ables Zi are mutually independent. The quantities it!j, (Tij and jump scaling factors Ji^i are con- 
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stant parameters. We write R = [-Ri, . . . ,Rn]', J = {■Ji,i)i=i,...,n;i=i,...,m, and cr = (o"ij)ij=i,...,n 
and we assume that the matrix S = crcr' is nonsingular. 

In Section [5], we will make further assumptions on the structure of the matrix S to facilitate 
the derivation of an explicit solution, assuming in particular that it possesses a factor structure. 
But the existence and structure of the optimal portfolio solution can be determined without 
further specialization, and we now turn to this problem. 



3.2 Wealth Dynamics and Expected Utility 



Let uo^t denote the percentage of wealth (or portfolio weight) invested at time t in the riskless 
asset and u!t = [(^i,t, ■ ■ ■ ,(^n,t]' denote the vector of portfolio weights in each of the n risky 
assets, assumed to be adapted caglad processes since the portfolio weights cannot anticipate 
the jumps. The portfolio weights satisfy 

En 
uji^t = 1- (3.4) 
1=1 

The investor consumes Ct at time t. In the absence of any income derived outside his 
investments in these assets, the investor's wealth, starting with the initial endowment Xq, 
follows the dynamics 



dXt = -Ctdt + uoiX, 



En 
1=1 



tXt 



dS, 



i,t 



= (rXt + cv'tRXt - Ct) dt + Xtu;[(TdWt + Xt J^^^ (0;^ J), Zi,tdNi,t. (3.5) 

We consider an investor with time-separable utility of consumption U{-) and subjective 
discount rate or "impatience" parameter /3 > 0. The investor's problem at any time t > is 
then to pick the consumption and portfolio weight processes {Cs, (ji>s}t<s<oo which maximize the 
infinite-horizon discounted expected utility of consumption. The optimal policies {Cs, <-t's}t<s<oo 
are subject to the admissibility condition that the discounted wealth process remains positive 
almost surely. 



Stochastic dynamic programming (see, e.g., [Fleming and Sonerl (120061 )) leads at time t to 
the discounted expected utility of consumption in the form V{Xt, Xt, t) where the value function 
is defined by 



V {x. A, t) = max Ea;,A,t 

{Cs,CJs; t<s<oo} 



-^'U{Cs)ds 



(3.6) 
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Here, the discounted wealth and intensities satisfy (13. 5 p and (13. 3 p over [t, oo) with initial con- 
ditions Xf = X, \t = A. Under the as sumption that V is sufficiently differentiable, the ap- 
propriate form of Ito's lemma (see, e.g., IprotteJ kooi )) for semi-martingale processes leads to 
the Hamilton- Jacobi-Bellman equation that characterizes the optimal solution to the investor's 
problem: 



= max 

{CM 



f dV{x,\,t) , ^ 

dt +Z."'(^ 
L 1=1 



- X,) + e-^^UiO 



Loo — 



(3.7) 



dV(x,\t), . ld'^V(x,X,t) 2 

ox 2 ox^ 

^ A/ J [Vix + {u)'J)i zx,\ + di,t)-V (x. A, t)] vi {dz) I 

with the transversality condition lim(_^oo E \V (X^, A(, t)] = 0. 

Using the standard time-homogeneity argument for infinite horizon problems, we have that 



e^^V {x, A, t) = max E^^x,t 

{C's,(^s; t<s<oo} 

= max E^^x,t 

{Cs,ijJs', t<S<Oo} 

= max Ex X o 

{Cs,u)s; 0<s<oo} ' ' 

= V^(x,A,0) = L(x, A) 



U{Cs)ds 



e-^^U{Ct^u)du 



e-^''U{Cu)du 



is independent of time. Thus V (x. A, t) = e ^^L{x, A) and (13. 7p reduces to the following time- 
independent equation for the value function L : 

= max I U{C) - /3L(x, A) + V a, (A,oo - A,) ^^^^^ 



+ 



dL (x. A) 

dx 



[rx + (jj'Rx — C) 



192L(x,A) 
2 



oj'T^ijjx^ 



m „ 

^ A/ / [L (x + {ijj'J)i zx, X + di) - L (x. A)] ui (dz) 



with the transversality condition 



(3.8) 



lim E [e-^'L (X^, A*)] = 0. 



t— >oo 



(3.9) 



12 



The maximization problem in (13 .Sp separates into one for C, with first order condition 

_ dL (x, A) 

^ - dx 

and one for uj : 



cj* = a;*(x, A) := argmax 



y ox 2 ax^ 

+ Ylk j [L{x + {uj'J)i zx, X + di)-L (x, A)] vi {dz) | (3.10) 

At time t > 0, given wealth and intensity vector Aj, the optimal consumption choice is 
therefore Q = C*{Xt, Xt) where 



C*(x, A) = [t/']"^ (5L (x, A) /9x) . (3.11) 

In order to determine the optimal portfolio weights, wealth and value function, we need to be 
more specific about the utility function U. 

4 Log-Utility Investors 

There are three classic utility functions for which one may hope to make further analytical 
progress, namely the log investor whose utility of consumption is the logarithm function, the 
power investor and the exponential investor. Collectively, these examples are known as HARA 
utilities. The optimal investment problem for various simpler types of market dynamics with 
these utilities lead to separable forms for the value functions. As we shall now show in this and 
Section [TJ this separation property extends to the Hawkes-diffusion model, albeit with some 
extra twists. In this section, we concentrate on the log-investor, for whom we are able to prove 
strong results on the existence and uniqueness of the optimal strategy. 

4.1 Optimal Investment with Log-Utility 

We now specialize the problem to that faced by an investor with logarithmic utility, U (x) = 
log (x). To start, we look for a candidate solution to (13. 8p in the form 

L(x,A) = /(A) + M-Mog(x) (4.1) 
13 



for some positive function / and constant M. Then 



^) _ f (X) ^) - M-V-i 



= -M-^x-l (4.2) 
and the optimal pohcy for the portfoho weights at time t > is o;^ = cj* (A^) where 
uj* (A) = argmin^ Ki{uj, A) 

Ki{uj,X) = l^-iv'R+^Lj'l^Lj-f^Xi j \og{l + {u'J),z)ui{dz)^. (4.3) 

We note that the convexity of Ki imphes this minimization has a unique solution u)* (A^) for 
any A G M". As for the optimal consumption policy, from equation (13. lip , and the facts that 
\U']~^{y) = y^^ and dL (x, X) /dx = M~^x~^, we obtain M = (3 and 

C; = (5Xt (4.4) 

Next, we substitute the optimal (C*,a;*) into (13.81) and determine that the function / must 
solve 

[^/](A)-/3/(A)=F(A), XeWl (4.5) 
where the Markov generator A for the process Af is given by 

m 

[•Am) = Yl (^'.- - W + A; [/ (A + d,)-f (A)]) (4.6) 



1=1 



and the nonhomogeneous term is 

F{X) = 1 - ^ - log/3 + {3-'Ki{u:*{X), A). 
The following lemma gives a computable formula for the smooth solution of (14. 5p . 

Lemma 1. The function f : — )■ M_|_ defined by 'is differentiable and given by the 

absolutely convergent integral 

POO 

/(A)= / e-^^Eo,A [F(A,)] rfs . (4.7) 
Jo 
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Proof. By the ergodic property, we know both that Eo,a [-^(-^s)] ^o,\ [-^('^oo)] and 

^Eo,A [-^(-^s)] — as s — > oo. One can also verify directly that there is a constant M > 

such that 

|F(A)| < M(l + ||Af ) . (4.8) 
From these facts follows the absolute convergence both of the integral in (14. 7p and the integral 

Of 



poo o 

/ e~^'—Eo,x[F{X,)]ds . 
Jo (J^t 



Since the right hand side of f l4.7p is differentiable, the Feynman-Kac formula implies it satisfies 
4.2 A Verification Result for the Log Investor 



he fo llowing verification theorem follows the logic outlined in Section III. 9 of 



Fleming and Soner 



(120061 ) and ensures that the above argument correctly characterizes both the optimal strategy 



and the associated value 
found in GoU and Kallsen 



unctio n. A more general verification result for log investors can be 

mm . 



Theorem 1. Consider the optimal problem Ii3.6\) for the log investor with impatience parameter 
(3 > 0, investing in the asset price model defined by liS. liS.S^) . liS. 3\) and satisfying the 
Stationarity Assumption that the matrix F = {ajSij — dij) is positive. 

1. The candidate solution V{x,\,t) = e"'^* [/(A) + (3~^\og{x)] is a classical (i.e. differen- 
tiable) solution of the HJB equation ([57 



2. For all initial conditions x > 0, A G M", the pair of processes (ci;*,C*),t > defined 
by ( [^.3[ ) and i4.4\) is an admissible policy, in the sense that they are progressively mea- 
surable and the process remains finite and positive {t, u) almost surely and solves the 
appropriate SDE. 

3. Let C denote the class of admissible policies {ujt, Ct),t > that satisfy 



limE,,A V{Xt,Xt,t) 

t—>-oc 



> 



(4.9) 
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For any {uj, C) G C, 



e-^'U{C,)ds 



< V{x,X, 0). 



(4.10) 



4- Let Vas denote the value function 

' POO 

Vas{x,\)= ^.,x / e-^-'UiC,)ds 

The optimal policy {u:l,C^ ),t > satisfies 



{A.U] 



lim E^^x 

t— >oo 



vix:,x„t) 



and the equality 



E 



e-^'U{Cl)ds 



y(a;,A,0). 



(4.12) 



(4.13) 



Hence V{x, A, 0) = Va5(x, A) and (a;*, C*) zs i/ie optimal portfolio policy in the class C. 



Proof. That V is a, classical solution of (13. 7p follows from Lemma [T] which shows that /(A) is a 
differentiable solution of ( 14. Sp . That {(jjl,C^),t > is admissible follows by general considera- 
tions. 

Suppose (cj, C) G C and consider the process = V{Xs, \s, s). For any t > 0, the Dynkin 
formula implies that 

E.,A[6]=eo+ f^.,x -(3V{X,, Xs, s) + f2»i (^^.oo - Am) 

J / 1 ^ 



dx 

m 

+ 5^A,,,_ 



+ ^ [rXs + uj^RXs - Cs) + a;,Sa;,X; (4.14) 



dx"^ 



V (X, + {uj'J)^ zX,, A,_ + di, s)-V (X„ A,_, s) ui (dz) 



ds 



Using the fact that V solves the HJB equation (13.71) leads in the usual way to the inequality 



E.,A [6] < Co - E,,A 



e-^'U{C,)ds 



(4.15) 



Finally, one can use the transversality condition ( 14. 9 p to take the limit t — > oo and obtain the 
desired result 



e~^'U{Cs)ds 



<^o = V{x,X,0). 
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The strategy {uj^,C*),t > maximizes (13 .Tp for almost every {t,uj), which means fl4.15p 
holds as an equality for every t > 0. Finally, one needs to verify (I4.12p in order to conclude 
that C,o = Vq = Vas- First, by the ergodic property fl2.5p one has 

hm E,,;, [e-^y{Xt)\] =0 

Next, by plugging in the optimal consumption C* = f3X* one then finds one can solve the 
problem of optimizing the expected utility of terminal wealth over the interval [0,t] with the 
adjusted interest rate f = r — (3 to show that 

E.,A [e-^*log(X;)] = [te-^*]4 fK,x[r + Ki{u;*{X,),Xs)]ds 

t Jo 

From the bounds (14.80 on the functions F, K this can be seen to converge to as t — )■ oo, again 
by the ergodic property (12. 5p . □ 

4.3 Properties of the Solution 

The log investor is often described as "myopic" because she acts at each moment in time as if 
the dynamical variables are in fact static parameters. Thus we should not be surprised that the 
optimal consumption and portfolio weights at any time given by (14. 4p and (14. 3 p are independent 
of the coefficients of the SDEs driving the asset returns. In particular, at any time t, the policy 
{ujI, C^, Xt) is precisely the same as the policy when the jump frequency is treated as a constant 
A = At, although that "constant" is changed at each instant. 

A second observation is that the policy (o;^, C^*, At), t > does not require knowledge of 
the function /(A). However, determining the function / requires solving the non-homogeneous 
equation (14. 5p . or equivalently evaluating the integral in (14. 7p . For this, we need some further 
structure on the problem, which we now add in order to derive the complete form of the solution. 

5 Additional Structure on the Diffusive and Jump Risks 

To find interesting examples of closed form optimal portfolio solutions, it is convenient to model 
the variance-covariance matrix E of the diffusive part of asset returns in such a way that its 
inverse is explicit. For this purpose, we adopt a modelling approach that is common in asset 
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pricing, namely to assume a factor structure. That is, we specify a block-structure for E 
consisting of k blocks of dimension m, with n — mk: 

/v,. ... \ 



nxn 



5^2,1 



-'2,m 



v 

with diagonal (or within-asset class) blocks 



-'m.m—l 



(5.1) 



^1,1 = vf 
kxk 



( 



\ 



1 Pl,l 
Pl,l ' ■ ■ 
• • • Pl,l 



Pl,l 
1 



(5.2) 



and off-diagonal (or across-asset class) blocks 



"^1,3 

kxk 







(5.3) 



where 1 > pi^i > — l/(/c — 1). 

The spectral decomposition of the E matrix is 



\— «,m 1 , v-^™ / 1 ,\ 



(5.4) 



where 

Ku^vf{l + {k-l)pi,i) (5.5) 
1^21 = vf (1 - pi^i) (5.6) 

are the 2m distinct eigenvalues of E. The multiplicity of each ku is 1, and the multiplicity of 
each K21 is k — 1. The eigenvector for kh is 1/, the n— vector with ones placed in the k rows 
corresponding to the /—block and zeros everywhere else, that is 

1^ = [0,...,0, 1^_^,0,...,0]', (5.7) 

asset class I 

where the first 1 is located in the k{l — + 1 coordinate. Mi is an n x n block diagonal matrix 
with a k X k identity matrix Ik placed in the Z— block and zeros everywhere else: 



Ml 

nxn 



••• 

; Ik ■ 
••• 



(5.8) 
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Corresponding to the above spectral structure, we have the orthogonal decomposition 
V © V-^ where V is the span of {1;};=!,..,™ and V-^ is the orthogonal space. 
As for the vector J of jump amplification coefficients, we assume that 



/ / T \ 

-'1,1 • • • 'Jl,m 



J — [Ji,...,J^ 

nxm 



Jn,m 



(5.9) 



where 



Ji ^ jih = [0^_^,...,jz^_^,...^0^^]' (5.10) 

asset class 1 asset class I asset class m 

for I — 1, ...,m. This structure means that assets within a given class I have the same response 
to the arrival of a jump, i.e., to a change in Nt- But the proportional response ji of assets of 
different classes to the arrival of a jump can be different {ji ^ jh for I ^ h). 
Finally, we assume that the vector of expected excess returns has the form 

Em _ . _ . 

^_^Rih + R = R + R^. (5.11) 

Here, we allow the expected excess returns to differ both within and across asset classes, by 
allowing R^ 7^ . The components of R play the role of the assets' alphas. The general R^ 
is orthogonal to each 1; and has the form 

R^ = [R-^ , R^\ 
where each of the /c- vectors R^ is orthogonal to the /c- vector 1. 

5.1 Closed- Form Optimal Portfolio Solution 

We now look for a vector of optimal portfolio weights w, and it is convenient to look for it in 
the form of its decomposition using the same basis as above. 



Em . . 

+ w^ = u> + w^ (5.12) 
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where u)^ = [ujf, ...,ujj"']'. The objective function K2{uj) in (I4.3P to be minimized reduces to 

1 1 1 



2^ ' "fc 

n ~ 

- Az / log (1 + kuiiiz) vi [dz) . 
1=1 ^ 

The minimization problem then separates as 

{u;^*,u:*) = argmin|^,^| {g^ (u,^) + g (u)} (5.13) 

where 

g^{u^) = -u^'R"^ + icu^'S^Lc;^ (5.14) 

(7(cj) = -/c^w;^; + -fc^wf/tiz - ^ Az,t / log(l + kuijiz) Vi (dz) . (5.15) 
1=1 1=1 1=1 

For the first part, minimizing g^{u:-^), the structure of S"*- implies that 

and therefore the optimal solution u)-^* has blocks 

o;,^* = —Ri (5.16) 

for / = 1, . . . , m. This part of the solution depends only on the diffusive characteristics (expected 
returns and variance-covariance) of the asset returns. 

The problem of minimizing g{(jj) separates itself into m separate minimization problems. 
With the change of variable 

Win = kui (5.17) 
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we see that 



S^-zuinRi + ^Tuf„Kii/k - Xi^t j log (1 + '(^injiz) ^1 {dz)^ (5.18) 



The convexity of the objective function imphes the existence of the minimizer. We can then 
determine in closed form. 

Two cases are exphcitly solvable. We first consider the case where the jump size is deter- 
ministic. In this situation, vi^dz) = 5 {z = zi) and the objective functions become: 

fn {W) = - 2^^^^ WlnRl + - 2^^^^ ^Inl^ll/k 
m 

- Xlt log (1 + ^InjlZl) . (5.19) 

1=1 

The first order conditions for the asset allocation parameters Ci7/„ are given by 

- Ri + winKii/k - Xi^tiiZi (1 + winiiZiY^ = for / = 1 , . . . , m. (5.20) 

These first order conditions form a system of m independent quadratic equations. Each separate 
equation f l5.20p admit a unique solution win satisfying the solvency constraint zuinjizi > —1. 
These are solvable in closed form: 



Wi -i^ii/k + jiZiRi + Y [jiZiRi + t^ii/k) + AXi^tjfzfnu/k 



A second case that is solvable in closed form is one where each jump term Zi^t has a binomial 
distribution {ui with probability pi or di with probability 1 — p;), since the corresponding first- 
order conditions are cubic. The first order conditions are obtained by differentiation with 
respect to ruin of the objective function stated in fl5.18p . namely: 

- ft + ^,„.„/* - A,,,„ / (1 + ^,„„.)-' z., (dr.) = for i = 1. . . . . m. (5.22) 

For binomially-distributed jumps, the conditions reduce to 

- Rl + WinKii/k - Xl^tjl {piui (1 + zuinjiui)'-^ + (1 " Pi) di (1 + TUinjldlY^) = (5.23) 

which produce a cubic polynomial equation in zuin, again exphcitly solvable, separately for each 
for / = 1, . . . , m. 
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6 Consequences for the Optimal Portfolio Allocation 



We now investigate in more detail the consequences of the exphcit portfoho weight formulae for 
an optimal asset allocation. The first element we note from f l5.2ip is the fact that the optimal 
portfolio is time- varying, since its composition changes with the jump intensities A^. 

In the case of purely diffusive risk, the optimal portfolio weights reduce to the classical 
Merton formula 



for / = 1, ..,m. As is well known, the solution in this case is constant. 

In the case where Poissonian jumps are added to the model, the solution specializes to fl5.2ip 
but with Xi^t replaced by the constant Poissonian jump intensity. As in the purely- diffusive case, 
the solution becomes constant. A solution with similar qualitative features would be obtained 
in the mutually exciting case if one replaced each stochastic jump intensity by its unconditional 
expected value, although given that the portfolio weights are nonlinear functions of Xi^t these 
would not be the unconditional expected values of the portfolio weights. 

Let us now return to the full solution in the mutually exciting case. A univariate model 
captures only part of the mutual excitation phenomenon: with a single asset, only time se- 
ries self-excitation can take place. In order to investigate the full potential impact of mutual 
excitation on optimal portfolio holdings, we now specialize the results above to a two-asset 
model where both time series and cross-sectional excitation can arise. Assets 1 and 2 can excite 
each other, not necessarily in a symmetric fashion, depending on the 2x2 matrix of mutually 
exciting intensities with coefficients dij, i,j = 1,2. The formulae above specialize with n = 2, 
k = 1 and m = 2, in which case Kh = vf for / = 1, 2 and hence: 



(6.1) 



or, replacing S by its explicit expression 




(6.2) 




(6.3) 
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Consider the change in the optimal portfoho allocation of an investor who observes a first 
shock, say to asset 1. For concreteness, let us return to the two-asset class scenario illustrated 
in Figure [3l In a Poissonian jump model, observing the first shock at time Ti does not change 
the investor's optimal portfolio: since jumps' future arrivals are independent of past jumps, 
there is nothing to do going forward other than to absorb the losses from the first jump. In 
the mutually exciting model, however, the occurrence of the first jump at time Ti self-excites 
the jump intensity Xi^t for t > Ti. This increase in Xi^t translates, if zi < 0, into a reduced 
asset allocation to asset 1. Moreover, these jumps have a contagious effect on ^2 since a jump 
in asset 1 cross-excites the jump intensity A2,t of asset 2. If Z2 < 0, then the optimal policy is 
to reduce the asset allocation to asset 2. Note that the reduction to the position in both risky 
assets occurs immediately after Ti, without waiting for future jumps. 

For the same sample paths as in Figure |31 Figure S] shows the optimal portfolio weights. 
This is a flight to quality, in the sense that the occurrence of a single jump in asset 1 causes 
the investor to flee both risky assets (or all of them in the general n case) for the safety of 
the riskless asset. This phenomenon is well documented as an empirical reality in practical 
situations; we believe that this is the first portfolio choice model to actually capture it in a 
theoretical setting. 

Increases in the jump intensities raise the probability of observing another jump in 5*1 at 
the future time T2. This, in turn, raises the probability of seeing a jump in S2 at time T3. 
Later on, at time T4, the jump in ^2 raises the probabihty of seeing a jump in Si at some 
future time T5, and so on. Mean reversion in the jump intensities at respective rates ai and 02 
counteracts these successive increases, keeping the intensities non-explosive (stationary, in fact). 
The optimal portfolio policy reacts to each change in jump intensity accordingly: increases lead 
to reduced asset allocation, decreases to increased asset allocation. Inevitably, this analysis is 
conducted in a partial equilibrium framework: it assumes among other things that expected 
returns do not change as jump intensities change, or at least not sufficiently to reverse the 
result. 

Another interesting empirical phenomenon that can be revisited in light of these optimal 
portfolio policies is home bias. Home bias refers to the observed tendency of most investors' 
portfolios to be insufficiently diversified internationally. In the model, the benefits from diversi- 
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Figure 4: Mutual excitation in a two asset-class world: Jump intensities (top panel) and optimal 
portfolio weights (bottom panel). 
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fication are much less valuable than in the standard diffusive model since international assets do 
not protect as much against jumps in domestic assets in the presence of cross-sectional mutual 
excitation. 

Finally, one phenomenon that is often documented in financial crises is the large increase 
in correlations between asset classes, with all of them increasing towards 1. In the context of 
the model, empirical correlations measured over a period where mutual excitation occurs will 
indeed be close to 1 as long as the jumps that result from mutual excitation are of the same 
sign (say both zi < and Z2 < 0), at least on average. In such periods, the jumps' contribution 
to the observed correlation trumps the continuous contribution. 

7 Other HARA Investors 

The cases of power and exponential utility also lead to candidate value functions and optimal 
portfolios in separable form. Due to the complexity of the underlying HJB conditions, the 
relevant verification result requires a lengthy, and perhaps uninformative, analysis that we have 
not yet completed. Nonetheless, as we now show, the candidate solutions can be characterized 
in terms of a fixed point problem that in principle can be solved numerically. 

7.1 Power Utility 

In this section, we consider the power investor with U{c) = c"'/'~f, 7 G (— oo,0) U (0,1). The 
analysis now consists in verifying the consistency of the following form for the solution to (13.81) 
in the form 

L(x, A) = x^g{\)/'y 
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for some positive function g. Substitution into (13. 8 p leads to 

= max {f/(C) - /3L (7.1^ 

{C,u}} 



J m r 
Ij 



ai {\oo - \i) + A; {g{X + di) - g{\)) 



+ 7L (r + u'R - C/x) + ^ ^''^ JHuj (7.2) 



+LY,hg{X + di)/g{X) / [(l + a;'j2)^-l]z/^(rf 



The optimal policy for the portfolio weight at time t > is a;* = ijj*{h{\t)) where 

J argmin^ir^(a;,A) 7>0 

a; (A) = <^ (7.3) 
I argmax^^ i^'^(LJ, A) 7<0 

hi{\) = \g{X + d{)/g{X) (7.4) 

K^(a;, A) = --iOj'R - ~ u'T.oj -f^^it [ [{I + uj' J zY - l] (dz) . (7.5) 

2 ' i(o,i] 



Solving the first order condition for C leads to the optimal consumption C* = x {'yg{X)) 
Finally, g is characterized by the implicit equation 



1/(7-1) 



[Ag]{X)-{f3-n)g{X) = G{X-g) (7.6) 
G(A; g) = g (A) K^' iu;*{h{X)), h{X)) + (1 - ^){^g{X)r/^''-'\ X G R!^. (7.7) 

The Markov generator A is again given by (14. 6p . We use the Feynman-Kac formula to write 
the solution g as a fixed point of an infinite dimensional nonlinear mapping: 

9 = Q{9) (7.8) 

poo 

{g{g)) (A) := / e-(^-^^)^Eo,A [G(A.; g)] ds. (7.9) 
Jo 

Remark 1. Note that the power case differs from the log case in two distinct ways. The 
first term on the right side of ( [7. 7| j is a distorted version of the right side of lji4-5\ ), where the 
A dependence in the function K is distorted in a g dependent fashion through the mapping 
h. The second term does not arise in the log case, and int roduces complications; a simila r 



situation occurs and is dealt with in a different model, see 
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Delona and Kliippelbergi (12008 ) 



Following xDelong and Kluppelhera ^200a ). one can attempt to verify that (TT^ is a contraction 
mapping, and that conseguently the seguence of iterates {g^^\i = 0,1,...} with g^'^^ = 1 and 
gi^+'i-) = converges to g. We do not attempt this here. 

In examples for which the functions K andu:l{X) are explicitly solvable, the iteration scheme 
can apparently be efficiently implemented numerically. 

7.2 Exponential Utility 

An investor with the exponential utihty U{x) = —e~^^j'y with risk aversion parameter 7 > 0, 
unhke the log investor, can in principle consume at a negative rate, perhaps even reaching 
negative wealth, and thus in this setting the question of defining admissible strategies is more 
involved. We can however, search for candidate optimal strategies that solve the reduced HJB 
equation f l3.8p in the form L(x, A) = — exp[— Ka;](y'(A) for some positive function (7, and then 
attempt to interpret the result. The HJB equation turns into: 

= max |f/(C) - /3L + -^\Ag\{\) 

{CM I g{\y 

- kL {rx + u;' Rx - C) + ^u;"SuJx^ (7.10) 

where Ag is given as before by ( 14. 6p . The first order conditions for C* imply that '~fU{C*) = kL 
and hence 

C* = — (kx — log g — log k) . 

7 

The candidate optimal portfolio weights are best expressed in terms of dollar amounts TTj^t = 
oJi^tXt invested. We find tt^ = 7r*(/i(A)), h = {hi, . . . , km) where 

7r*(A) = argmin^ K{n, A)) (7.11) 
hi{X) = Xig{X + di)/g{X) (7.12) 

2 in „ 

K{tz, A) = KTz'R - yTr'STT -^>^i (exp[-/t (tt' J), z] - 1) vi [dz) 

1=1 

Substitution of C*,7r* back into f l7.10p leads to 

= (r - /3)^ + [^^](A) -K{rx- C*) - g{X)K{7r*{X), h{X)) 
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which in turn imphes that k — r^. The condition on g is now imphcit: 



[Ag]{X) -iP-n + rlog{n))g{\) = ^(A; g) 



(7.13) 



G(A; g) = g{X) [r logg + K{n*{X), h{X))] 



(7.14) 



This is very similar to the characterization of g for the power utihty case. In examples where 
K and tt* are explicitly known, one can attempt to solve numerically for g by iteration. 



This paper extends the range of models for which solutions to the optimal dynamic portfolio- 
consumption problem are available, to one which includes mutually exciting jumps. We analyze 
features of the optimal solution and show that it differs from the usual case in important ways. 
In particular, it introduces an explicit time- variation in the optimal portfolio weights in response 
to changes in the jump intensity, providing a rare example of an explicit time- varying optimal 
portfolio solution. Moreover, power and exponential investors both adopt more aggressive 
strategies, characterized by a distortion function h, than the corresponding investor who does 
not fully recognize the mutual excitation effect. 



8 Conclusions 
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